Measuring the Nematic Order of Colloidal fd Virus by X-ray Diffraction 
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The orientational distribution function of the nematic phase of the semi-flexible rod-like virus fd 
is measured by x-ray diffraction as a function of concentration and ionic strength. X-ray diffraction 
from a single-domain nematic phase of fd is influenced by interparticle correlations at low angle 
while only intraparticle scatter contributes at high angle. Consequently, the angular distribution 
of the scattered intensity arises from only the single particle orientational distribution function 
at high angle but it also includes spatial and orientational correlations at low angle. Experimental 
measurements of the orientational distribution function from both the interparticle and intraparticle 
scattering were made to test whether the correlations present in interparticle scatter influence the 
measurement of the single particle orientational distribution function. It was found that the two 
types of scatter yield consistent values for the nematic order parameter. It was also found that 
x-ray diffraction is insensitive to the orientational distribution function's precise form, and the 
measured angular intensity distribution is described equally well by both Onsager's trial function 
and a Gaussian. At high ionic strength the order parameter S of the nematic phase coexisting 
with the isotropic phase approaches theoretical predictions for long semi-flexible rods S = 0.55, but 
deviations from theory increase with decreasing ionic strength. The concentration dependence of 
the nematic order parameter was also found to better agree with theoretical predictions at high ionic 
strength, indicating that electrostatic interactions have a measurable effect on the nematic order 
parameter. The measured x-ray order parameters are also shown to be proportional to the measured 
birefringence and the saturation birefringence of fd is measured, enabling a simple, inexpensive way 
to measure the order parameter. 



I. INTRODUCTION 

The first step to understand the effect of steric inter- 
actions in colloidal rod systems was taken in 1949, when 
Onsager wrote his seminal paper on the phase behavior 
of hard and charged rods jlj . Onsager developed a free 
energy theory at the second virial level describing the 
phase transition of hard rods from an isotropic phase, in 
which the particles are randomly oriented, to a nematic 
phase, in which particles are oriented in a distribution 
about a preferred direction. All theoretical predictions 
for the properties of this phase transition, such as the 
coexistence concentrations and the nematic order param- 
eter, depend on the functional form of the orientational 
distribution of the rods in the nematic phase. Onsager 
chose one test function and in a later review paper Odijk 
showed that qualitatively similar results for the proper- 
ties of the phase transition can be found by choosing a 
Gaussian test function [^j. The exact form of the orien- 
tational distribution function that satisfies the Onsager 
theory can be obtained via series expansion j| [| [| or 
by direct iterative methods j^, f?|. Determining the ori- 
entational distribution function of the nematic phase of a 
colloidal rod system is the most sensitive test of whether 
a system is described by Onsager's theory. 

In this paper we measure the measure the concentra- 
tion and ionic strength dependence of the orientational 
distribution function of fd virus via x-ray diffraction. The 
fd virus is a charged semi-flexible rod with a length L to 



diameter D ratio L/D s» 130. The electrostatic charge on 
the rods can be taken into account by an effective diame- 
ter D e ff, larger than the bare diameter, which is approxi- 
mately equal to the distance between particles when the 
interaction potential is about k^T . The exact calculation 
for the effective diameter is outlined in Ref. [|, |J ■ An 
increase in ionic strength of the solution containing the 
charged rods produces a decrease in effective diameter. 
In Onsager's theory, the limit of stability of the isotropic 
phase is predicted to be ir/AnD c sL 2 — 4, where n is 
the number density Q. This is predicted to be valid for 
long rods with a length to effective diameter ratio greater 
than 100 p0| . Previously the isotropic and nematic coex- 
istence concentrations fd have been measured as a func- 
tion of ionic strength and it has been shown to agree well 
with numerical results from Chen for a semi-flexible rod 
of with a ratio of persistence to contour length of 2.5 
0. Theoretical models suggest that semi- flexibility also 
acts to significantly lower the nematic order parameter 
at coexistence. For fd, a relatively rigid polymer with 
a ratio of persistence to contour length of 2.5, the ne- 
matic order parameter at coexistence is predicted to be 
S = 0.55, which is significantly smaller than predicted 
for rigid rods, approximately S — 0.79 pT[ . Several re- 
view articles describe in more detail the theoretical and 
experimental aspects of this and other systems described 
by Onsager's theory 

In x-ray diffraction, the scattered intensity consists 
of two parts, intraparticle scatter F(q) and interparticle 
scatter S(q). The intensity can be written as a product 
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of the two types of scatter 



I(q)=NF(q)S(qJ(q)) 



(1) 



where q is the three dimensional reciprocal vector in 
cylindrical coordinates q = (q r , q z , <ft). In a uniaxial ne- 
matic, q r is perpendicular to the nematic director and 
the scattered intensity is independent of the azimuthal 
angle <f> about the director. If the system is oriented such 
that the nematic director is in the z direction, q can be 
described by q = (q r ,q z ). The intraparticle interference, 
or form factor, contains information about the structure 
of the individual particles. F(q) can also be written as 
< f(o) 2 > where f(q) is the fourier transform of the 
electron density of a particle and the average is over all 
the particles and their orientations. The interparticle 
interference, or structure factor term, contains informa- 
tion about the positional and orientational correlations 
between particles. The structure factor depends on the 
positions of the centers of gravity of two scatterers Ri, Rj 
and their relative orientations |l9fl: 
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The orientation of the particles is included in f(q) and 
the average <> is over all particles and their orienta- 
tions. For isotropic scatterers, fi(q) = fj(q) and the 
structure factor and the form factor decouple, but for 
anisotropic scatterers, fi(q) ^ fj(q) unless the particle 
orientations are the same. Therefore, in contrast to scat- 
ter from spheres, the structure factor S(q) of rods can 
not, in general, be decoupled from its anisotropic form 
factor F(q). 

In a nematic system, however, there is no long ranged 
translational order. As a result, S(q) approaches unity in 
the limit of high q, and if S(q) = 1, the scattered intensity 
is due only to the intraparticle interference diffraction 
and I(q) = F{q). In this regime the angular distribution 
of the scattered intensity is a function only of the single 
particle orientational distribution function. Because of 
the crystalline internal structure of viruses such as fd and 
Tobacco Mosaic Virus (TMV) , x-ray diffraction produces 
a complex pattern of intraparticle scatter at high q which 
can be used to measure the single particle orientational 
distribution function of the viruses p0| . 

At low q the scattered intensity is dominated by S(q), 
and the angular distribution of the interparticle inter- 
ference scatter is influenced by the angular and spatial 
correlations between neighboring rods. When intraparti- 
cle interference scattering is absent or too weak to inter- 
pret, as in thermotropic liquid crystal systems pil , or 
the system of lyotropic vanadium pentoxide (V2O5) Jl8| , 
x-ray investigations of the nematic orientational distri- 
bution rely on measuring the angular distribution from 
interparticle interference scattering. In this case one does 
not calculate the single particle orientational distribution 
function, but instead the coupled fluctuations of neigh- 



boring rods; this is predicted to overestimate the value of 
the nematic order parameter for highly ordered samples 

In this paper we explore the behavior of the nematic 
phase of fd virus, investigating the concentration and 
ionic strength dependance of the spatial and orienta- 
tional ordering. We present measurements of the orienta- 
tional ordering of the nematic phase in coexistence with 
the isotropic phase as a function of ionic strength and 
compare the results with the predictions for semi-flexible 
rods. Previously, measurements of the orientational dis- 
tribution function of a nematic phase have been made 
either from form factor scatter as in work done by Old- 
enbourg et al on TMV Q and work done by Groot et 
al and Kassapidou et al on persistence lengthed DNA 
fragments ]l7[ pi[ or from structure factor scatter as in 
work done by Davidson et al [|l8| . Using fd as our model 
rod allows us to not only measure the orientational dis- 
tribution function from intraparticle scattering, but also 
from interparticle interference scattering. This allows us 
to experimentally resolve this question of whether or not 
correlations between angular and spatial order present 
in interparticle scatter influence the measurement of the 
order parameter. The saturation birefringence of fd was 
also measured allowing for measurements of the order pa- 
rameter to be made using birefringence methods which 
involve much simpler techniques than x-ray diffraction. 

This paper is organized in the following manner. In 
section |0| we describe the virus system and the experi- 
mental methods. In section III qualitative observations 



about the diffraction data are made. This is followed by 
a description of the analysis technique used to extract 
the orientational distribution function from the diffrac- 
tion data in section |IV[ . Quantitative measurements of 
the nematic spatial ordering and orientational ordering 
are presented in section [y| This includes first a brief sec- 
tion describing the measured spatial ordering and then 
a section presenting the measured orientational ordering 
of the nematic fd. Section [v| summarizes the significant 
results of this paper. 



II. MATERIALS AND METHODS 

The physical characteristics of the bacteriophage fd are 
its length L— 880 nm, diameter D— 6.6 nm, persistence 
length p= 2200 nm and charge per unit length of around 
I0e - / nm at pH 8.2 [|l5|. When in solution, fd exhibits 
isotropic, cholesteric, and smectic phases with increasing 
concentration [^5], |6[ |2?|, ^] . The fd virus was prepared 
using standard biological protocols found in Ref. [g9| us- 
ing the JM10I strain of E. coli as the host bacteria. The 
standard yield is «50 mg of fd per liter of infected bac- 
teria, and virus is typically grown in 10-12 liter batches. 
The virus was extensively dialyzed against a 20 mM Tris- 
HC1 buffer at pH 8.2 and the ionic strength was adjusted 
by adding NaCl. 

X-ray diffraction was done at the SAXS station on 
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beamline 8-ID at the Advanced Photon Source at Ar- 
gonne National Lab. The beam flux is 2 x 10 10 pho- 
tons/s for a 50 x 50^m beam with a photon energy of 
7.664 KeV (A=1.617 A). The samples were a suspen- 
sion of monodisperse fd in the cholesteric phase, sealed 
in « 0.7 mm quartz x-ray capillaries. Cholesteric sam- 
ples were unwound and aligned in a 2 T permanent mag- 
net (SAM-2 Hummingbird Instruments, Arlington, MA 
02474)^3), forming a single domain nematic phase par- 
allel to the long axis of the capillary and the magnetic 
field, which we will call z. The free energy difference 
between the cholesteric and nematic phases is negligi- 
ble, and the theory of the phase behavior of the isotropic 
to nematic transition can be applied equally well to the 
isotropic to cholesteric transition observed in fd || . The 
magnetic field does not have a significant effect on the 
ordering of the nematic phase ||l], [32| . Samples had 
to remain in the magnetic field for a minimum of 15 min- 
utes at low concentrations and a maximum of about 8 
hours at the highest concentrations. The strength of 
the magnetic field limited the maximum concentration 
at which we could unwind the cholesteric phase into a 
mono-domain nematic to about 100 mg/ml [^]. Align- 
ment of the nematic sample was checked in a polarizing 
microscope and, using a 3A Berek compensator, its op- 
tical retardance was measured. To easily view the solu- 
tion within the capillaries under the microscope, samples 
were placed in an index matching water bath while still 
within the magnet. Birefringence is related to optical re- 
tardance R by An = R/d; d is the measured thickness 
of the nematic sample within the capillaries. The mag- 
net and sample were then mounted in a vacuum chamber 
such that the sample was in the beam line, and the mag- 
netic field was perpendicular to the incoming beam. To 
observe the effect of charge on the nematic phase, sam- 
ples were prepared at different concentrations and ionic 
strengths. The fd concentration was measured with a UV 
spectrometer by absorption at 269 nm with an absorption 
coefficient of 3.84 cm 2 mg _1 . 

When the solutions of fd were exposed to x-rays for 
extended time, disclination lines that matched the pat- 
tern traced by the beam could be seen with a polarizing 
microscope. Since our samples were exposed for varied 
times, a series of x-ray diffraction patterns from the sam- 
ples were collected with increasing x-ray exposure time 
to quantify sample damage and its effects on the scatter- 
ing pattern. The polarizing microscope revealed sample 
changes after w 6 s of exposure, but the angular spread 
of the diffraction peaks was not affected until exposure 
times increased above 10 s, at which point the angular 
interference peak scatter broadened significantly. The ef- 
fect of exposure for < 10 s on the calculation of the order 
parameter was not measurable. Data was collected for 
the interparticle interference scatter by averaging ten 10 
s exposures taken at different 50 x 50um sections. To 
observe the much less intense intraparticle scatter, the 
sample was continuously moved through the 50 x 50um 
beam allowing for a total exposure of 120 s. A single long 



exposure was used to image intraparticle scatter because 
it resulted in less noise than multiple short exposures 
because readout noise on the CCD was higher than the 
dark current. Readout noise and solvent scatter were 
subtracted from data images during analysis, but over 
the q-range which was analyzed this background scatter 
was very uniform and could be approximated as a con- 
stant. 



III. OBSERVATIONS 

The two dimensional scattered intensity of low an- 
gle interparticle and high angle intraparticle interference 
peaks are shown in Fig. ^ for concentrations spanning the 
range over which fd is nematic at 10 mM ionic strength. 
The angular spread of both types of scatter broadens with 
decreasing fd concentration or increasing ionic strength, 
corresponding to an increase in disorientation of the rods. 
The low angle peak, in Fig. |l|a, is very intense and is a 
result of interparticle interference. The horizontal posi- 
tion of this peak is inversely proportional to the average 
interparticle separation, and the radial width of the peak 
is inversely proportional to the correlation length of the 
interacting rods. At larger scattering angle, the fd layer 
lines are visible as shown in Fig. [l]b. These intraparticle 
peaks are much less intense than the interparticle inter- 
ference peaks and are the result of single particle scatter 
arising from the helical packing of the virial proteins. 
The layer lines occur at intervals along the z direction 
proportional to the reciprocal of the axial repeat of the 
helical protein coat, which is 33Aj3j| . Because of discrep- 
ancy in both intensities and scattering angle between the 
interparticle and intraparticle scatter, we were unable to 
image both the high and low angle scatter simultane- 
ously. 

Because of the short ranged positional order in the ne- 
matic phase these intraparticle interference peaks should 
be independent of interparticle correlations. We con- 
firmed this hypothesis by observing that the location of 
the peaks do not change with concentration as the inter- 
particle peaks do. We also compared our data to pub- 
lished fiber diffraction results for M13[§|, Ml 3 is 
also a filamentous bacteriophage, which only differs from 
fd by one amino acid per coat protein: their structures 
are otherwise identical and virtually indistinguishable by 
x-ray fiber diffraction |56j. Upon comparing published 
fiber diffraction data with our data from nematic fd, 
we observed that they were similar, but that the fiber 
diffraction patterns had Bragg peaks due to the hexag- 
onal packing of the virus in the fiber which were absent 
in our nematic diffraction data. We also noticed that the 
horizontal location of the single particle peaks in the fiber 
diffraction was 4% larger than the location of our solu- 
tion diffraction peaks, indicating that the fiber diffraction 
was done on virus which had a smaller diameter than 
those in our nematic samples. The fibers are partially 
dehydrated, so it is not surprising that they become com- 
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FIG. 1: (a) Contour plots of scattering from nematic fd due 
to interparticle scatter, (b) Contour plots of scattering from 
nematic fd due to intraparticle interference. The interparticle 
scatter shown in (a) is hidden behind the beamstop in (b) 
which is located on the left of the images and can not be seen 
in the contour plot. From top to bottom the concentration 
of the samples are 93 mg/ml, 33 mg/ml, and 15.5 mg/ml. 
Samples shown are at an ionic strength of 10 mM (20mM 
Tris buffer), pH 8.2. The magnetic field and virus orientation 
are perpendicular to the scatter as shown in the schematic. 
$ is the angle from the equator on the detector film. Axes of 
contour plot are labeled in A -1 , note the scales are different 
in (a) and (b). 



pressed. The layer line spacing, however, was not altered, 
indicating that no stretching of the virus occurs in the 
fibers. From these observations we concluded that the 
high angle scatter from the nematic fd was independent 
of interparticle correlations. Detailed analysis of both the 
interparticle and intraparticle diffraction continues in the 
following sections. 



IV. ORIENTATIONAL ANALYSIS TECHNIQUE 



Because of the short ranged positional order of the 
nematic phase, the high angle scattered intensity should 
be independent of interparticle correlations S(q r ,q z ) = 1. 
We have demonstrated above that this is true for fd. fn 
this case, the intraparticle scattered intensity of a system 
of rods is related to the orientation of those rods in the 



following manner [£0|, p7[ : 

I(q r ,q z ) =< I s (q r ,q z ) >= J $(0)/ s (? r [0], g,[fi])dn (3) 

Where is the solid angle (8, <fr) a rod makes with respect 
to the nematic director 9 and azimuthally with respect 
to the incident beam <j>. Because fd is axially symmetric 
$(57) simplifies to $(0). I s (q r , Qz) 1S the axially symmet- 
ric three dimensional form factor f(q) 2 of a single rod. 
<&(#) is the orientational distribution function (ODF) of 
the rods. Because the form of the ODF is not known 
exactly, three test functions were used: 



$(0)= Aexp-^ (0<e<n/2) 
= Aexp-^g^ (%/2<6<%) , 
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where a sets the width of each of the peaked func- 
tions, and A is the normalization constant such that 
/ $(0)sin(0)d0# = 1. Eq. | is normalized. The first 
ODF is the Gaussian used by Odijk [EL the second is the 
function used by Oldenbourg et al. P0| in their study of 
diffraction from nematic TMV, and the third was defined 
by Onsager. The nematic order parameter, 
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£ Qcos 2 (0)-i)$(tf)dcos(#) 



(7) 



was determined for the orientational distribution func- 
tions which best described the diffraction patterns. 

The scatter from intraparticle interference was ana- 
lyzed by comparing it to a simulated scatter created from 
the evaluation of Eq. ^| using a three dimensional model 
for the single rod form factor and a trial ODF. Previ- 
ously Oldenbourg et. al. measured the ODF from the 
intraparticle interference scatter of TMV by simplifying 
Eq. H to a one dimensional integral at a constant q r [|o) . 
This one dimensional method could not be used for in- 
traparticle fd scatter because fd has a protein coat with 
a pitch much larger than that of TMV, 33A versus 23A 
respectively, resulting in layer line overlap at low con- 
centrations. Instead, the radial intensity distribution of 
single rod was modeled by B7| 



(8) 



I s (q r ,q z ) = I m {q r ,q z )W2'Kaq r 



I m is the scattered intensities along the middle of the ze- 
roth and ± first layer lines of our most aligned nematic 
sample, S — 0.96 and Gaussian a — 0.11 as determined 
by the interparticle interference peak. The intraparticle 
interference data that fell on the detector in the range 
of q r = 0.19 - 0.33 A, which encompasses the innermost 
peak on each of the three layer lines visible in the interfer- 
ence pattern, was fitted to the model diffraction images. 
For each diffraction pattern, an a was found for each 
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trial distribution function which minimized a computed 
chi-squared value 
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where B and C are fitting parameters and i sums over the 
pixels in the scattered image. B was calculated once for 
each scattered image, and was not adjusted when com- 
paring different ODF's. For more details of the model 
and analysis of the intraparticle diffraction images, refer 
to Appendix [a|. 

To measure the orientational distribution function 
from the interparticle peak, the method of Oldenbourg 
et al. was used, because the scatter consists of only one 
peak. In this method, Equation [| simplifies to a one di- 
mensional integral at constant q,.. This method is identi- 
cal to that frequently used for analyzing thermotropic in- 
terparticle scatter, with the exception that Oldenbourg's 
method approximates the single rod scattering, I s (f2) , as 
being proportional to l/sin(w) for small 9. In this equa- 
tion, lj is the angle between the rod and the incident 
beam. This is in contrast to most analysis done on in- 
terparticle interference in the past, in which / s (^) = 1 
pjl |22| , The l/sin(w) proportionality attempts to 
include finite size of the rod into the calculation of the 
ODF. For more details refer to Appendix |X| 



V. RESULTS AND DISCUSSION 

A. Spatial Ordering 

The location of the maximum, q m , of the first inter- 
ference peak and its radial width, Aq m , were measured 
along the equator, q z = in order to obtain information 
about the spatial ordering of the system. Because we are 
only analyzing data along the equator, these properties 
can be determined by dividing the equatorial form fac- 
tor, F(q r , 0), from the scattered intensity peaks, l(q r ,0), 
and then by fitting the remaining structure factor peak, 
S(q r ,0), to a Saussian S(q r ,0) = e -(9™-?.) 2 /2(A? m ) 2 as 
done in Ref. |3!| . I(q r , 0) and S(q r , 0) are shown in Fig. 
U for three different samples. The equatorial form fac- 
tor scatter was approximated by the Fourier transform 
of the known equatorial projection of the cylindrically 
averaged electron density of fd |Q . The electron density 
was approximated by binning the radial electron density 
into 10 sections as illustrated in the inset of Fig. ||a. The 
location of the equatorial peaks produced by the Fourier 
transform of the electron density agree with the equato- 
rial form factor data obtained at higher angle, but the 
increase in S(q r ,0) at high q r shown in Fig. |^ indicates 
that this approximation is only qualitatively correct at 
high q r and that the presence of background noise in the 
interparticle diffraction data hides any high q r form fac- 
tor information. At high concentrations the scattered 
intensity is much stronger compared to the readout noise 
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FIG. 2: Equatorial intensity profile, I(q r ), and equatorial 
structure factor, S(q r ), for three representative samples at 
10 mM ionic strength pH 8.2. Smaller inset graph is the 
binned radially averaged electron density p used to calculate 
the equatorial form factor shown as a dashed line. The de- 
viation of the structure factor from one at high q r is due to 
both background noise in I(q r ), which hides the actual form 
factor and a loss of accuracy in the model form factor at high 
gr- 



and as a result we are able to analyze the structure factor 
data to higher q r than at low concentrations. 

The q m and Aq m measured are plotted as a function 
of concentration for two different ionic strengths in Fig. 
||a. With increasing concentration, the average rod sep- 
aration decreases as c -1 / 2 (q m oc c 1 ^ 2 ) as expected for 
any rod system pS| . At a given concentration the rod 
separation remains constant and the variance increases 
with decreasing ionic strength. The electrostatic repul- 
sion present between the rods causes the rods to main- 
tain the maximum separation possible, but a smaller ef- 
fective diameter at low ionic strength allows for more 
fluctuations. The number of rods per correlation length 
q m /Aq m is plotted as a function of concentration in Fig. 
|3|b. The concentration dependance of q m j Aq m is much 
more significant at 10 mM ionic strength, than at 110 
mM, indicating that at high ionic strength the rods are 
less correlated. It is interesting to note that the second 
interference peak is much weaker than the first interfer- 
ence peak, indicating a large Debye- Waller factor. This 
is in contrast to charged 3D spherical and 2D disk sys- 
tems which show a much stronger second, and even third 
interference peak |ll], The structure factor of nc- 

matic fd also contrasts that of nematic end-to-end ag- 
gregated TMV, a very rigid rod, which has a structure 
factor closely resembling that of the 2D disk systems [[39| . 
One way to interpret the large and sharp first peak in the 
structure factor of fd is that flexible nematic rods have 
long range spatial correlations similar to a dense fluid of 
disks. However the absence of secondary peaks in the 
structure factor implies that fd particles have a greater 
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FIG. 3: (a) The concentration dependance of the maximum 
q m of the interparticle interference peak. The average rod 
separation is a distance of 2-7r/q m A. The equation of the curve 
fitted to the combined data sets is q m — 0.004c 1//2 . The inset 
graph shows concentration dependance of the variance of the 
interference peaks Aq m . (b) The concentration dependance of 
q m / Aq m the number of rods per correlation length. Squares 
(□) are at 10 mM and triangles (T) are at 110 mM ionic 
strength pH 8.2. 



degree of positional disorder about their average position 
than do disks. Perhaps the flexibility of fd accounts for 
this dramatic difference in spatial organization. 



B. Orientational Ordering 

By examining the \ 2 values obtained from the inter- 
and intra-particle scatter analysis, and the residues 
(Idata-I/it) from the interparticle scatter analysis we de- 
termined that analysis of x-ray diffraction data does not 
yield a unique orientational distribution function. The 
Gaussian and the Onsager distribution function each fit 
the intensity data equally well when comparing residues 
and x 2 values from each of the two functions. However, 
we were able to eliminate Oldenbourg's distribution func- 
tion from the possible ODF forms because it did not ac- 
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FIG. 4: (a) Angular intensity scan at g r =0.07 ± O.OOlA -1 
from the three data scatter shown in Fig. fi]a with best-fit 
curves calculated from the three trial ODF. Solid fit line rep- 
resents fit of both the Gaussian and Onsager ODF's, dotted 
line is the fit of Oldenbourg's ODF. (b) Residue (Idata -I fa) 
plot. 9 is illustrated in Fig. Fla. 



curately model the tails of the diffraction data at low 
concentration. This insensitivity of x-ray diffraction to 
the exact form of the ODF was predicted by Hamley 
who showed that x-ray patterns are insensitive to higher 
order terms in the spherical harmonic expansion of the 
orientational distribution function and therefore only an 
approximation to the full orientational distribution func- 
tion can be found Q . 

To demonstrate this assertion, the scattered interpar- 
ticle intensity at a constant radius of q r = 0.07 ± 0.001 
A -1 is plotted in Fig. [|a with the best-fit model intensi- 
ties for each of the three ODFs. 'J is the angle from the 
equator on the detector film. The actual best-fit orienta- 
tional distribution functions calculated from these inter- 
particle angular scans are shown in Fig. ||. The residues 
calculated from the interparticle and intraparticle inter- 
ference results for the three samples are illustrated in 
Figs. |^b and ||c, respectively. The intraparticle scatter 
residues shown are for the scattered intensity shown in ^a 
minus the model images shown in Fig. created with 
the Gaussian ODF. The intraparticle model scatter pro- 
duced relatively uniform residues indicating that it was 
a qualitatively good model. In two dimensions, we were 
unable to distinguish differences between residue plots of 
ODFs of the same width, therefore residue analysis was 
limited to the interparticle scatter. At high concentra- 
tion small systematic disagreements between the best-fit 
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FIG. 5: Orientational distribution functions calculated from 
the interparticle angular intensity scan at constant radius 
q r =0.07 ± O.OOlA -1 shown in Fig. §b. Gaussian (thin solid 
line) , Oldenbourg (dotted line) and Onsager (thick solid line) 
ODF are shown. Order parameters shown are calculated from 
each ODF. From top to bottom the concentration of the sam- 
ples are 93 mg/ml, 33 mg/ml, and 15.5 mg/ml. The ionic 
strength of the samples is 10 mM, pH 8.2. Axes of contour 
plot are labeled in A" 1 . 



models and the data are most visible in the residue plots 
in Figs. and ^c, but each of the three models and 
their respective ODFs are nearly indistinguishable. Ex- 
cept at low concentration, the best-fit model intensities 
obtained from the three distribution functions can not be 
distinguished from one another both by analyzing residue 
plots and by comparing minimum \ 2 values computed 
from the fitting routine. At low concentration the sys- 
tematic disagreements between the data and the fits are 
lost in the noise, but disagreements in fits from differ- 
ent ODFs become visible. The best fit model intensities 
from the Gaussian and Onsager ODFs are indistinguish- 
able, but the residues from the Oldenbourg ODF show 
disagreement, and the fits are systematically higher than 




FIG. 6: (a) Contour plots from Fig. |ljt> of scattering from 
nematic fd samples due to intraparticle interference. Bottom 
scatter is at 15.5 mg/ml and top is at 93 mg/ml. Samples were 
at 10 mM ionic strength pH 8.2. (b) Simulated intraparticle 
scatter using a Gaussian ODF which best fit the intraparti- 
cle scatter shown in (a), (c) Residue (Idata — 1 fit) /l fit plot. 
Maximum residues in (c) are ±10%. The axis of the plots are 
labeled in A -1 . 



the background scatter at high angle '5 . At the isotropic- 
nematic transition the \ 2 values computed from the Old- 
enbourg ODF were also consistently higher. The calcu- 
lated Oldenbourg ODF also looks significantly different 
from the calculated Gaussian and Onsager ODFs. From 
these qualitative observations we argue that the distri- 
bution function used by Oldenbourg et al. does not de- 
scribe our diffraction data as well as the Gaussian and 
the Onsager distribution function at low concentrations. 
The Gaussian and the Onsager orientational distribution 
functions fit the diffraction data equally well. 

Because of small differences in the trial orientational 
distribution functions (as illustrated in Fig. |J), best-fit 
ODFs vary slightly in their width, and subsequently re- 
turned slightly different order parameters. But, the order 
parameters calculated from the best-fit Gaussian and On- 
sager ODFs were in agreement with one another within 
the experimental uncertainty of AS/ S < 6%. Order pa- 
rameters calculated with the Oldenbourg ODF were in 
common agreement at high concentrations, where model 
scatter agreed with the data. The nematic order param- 
eter calculated at multiple q r across interparticle peak 
also remained relatively constant, AS/S < 4%. Because 
we can not distinguish between the Gaussian and the 
Onsager model scatter, the order parameters to be pre- 
sented henceforth are an average of the values calculated 
from only the Gaussian and the Onsager ODF, and the 
uncertainty on the values given are a combination of ex- 
perimental error and uncertainty due to variation in order 
parameters from two trial ODFs. 

The concentration dependance of the order parameters 
was measured from both the interparticle and intraparti- 
cle peaks and the resulting values are graphed in Fig. 
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[?]. In Fig. U the order parameter of the nematic phase 
in coexistence with the isotropic phase is plotted for five 
different ionic strengths as a function of concentration. 
The coexistence concentrations are an increasing func- 
tion of ionic strength. Our analysis shows that the order 
parameters calculated from the interparticle and intra- 
particle scatter are consistent with one another both as 
a function of concentration and of ionic strength, indi- 
cating that correlations in the interparticle peak do not 
visibly change measured nematic order parameters. 

Fig. |^a shows data obtained at an ionic strength of 
10 mM and Fig. ^|b shows data obtained at an ionic 
strength of 110 mM, pH 8.2. With increasing concen- 
tration, the order parameter increases until it saturates 
near S = 1, and at constant concentration, the nematic 
order parameter decreases with increasing ionic strength. 
At low concentrations, the scattered intensity is spread 
over a large area due to the broad orientational distri- 
bution function, which leads to a large decrease in the 
signal to noise ratio. As a result, the variation in the 
calculated order parameters increases, reaching a maxi- 
mum of AS/S < 10%. The theoretical curves shown in 
Fig. ^ were computed from a scaled-particle theory which 
includes semi-flexibility in the orientational entropy and 
electrostatic interactions through an effective diameter. 
This calculation is outlined in detail in Appendix [§]. The 
Onsager ODF was used in calculating these theoretical 
curves. Our results qualitatively agree with this theory 
at low ionic strength, and quantitatively agree at high 
ionic strength. 

The values for the order parameter calculated from the 
x-ray diffraction data were also compared to birefringence 
measurements for the whole range of concentrations and 
the results are plotted in Fig. u. Theoretically we expect 



An _ Ansat 
c c 



(10) 



where An is the sample birefringence and An sat is the 
birefringence of perfectly aligned fd [Q. The satura- 
tion birefringence per unit concentration was measured as 
An sat /c = 3.8 x 10~ 5 ±0.3 x 10~ 5 ml/mg using data from 
samples at five different ionic strengths. Our value for the 
saturation birefringence is lower than the previously cal- 
culated value by Torbet et. al. of An sat /c = 6 x 10~ 5 
ml/mg which was calculated for solutions of fd at 16 
mg/ml in 10 mM Tris-HCL buffer at pH 7.5 placed in 
a 2-4 T magnetic field @. 

For long rods Onsager predicts that the nematic order 
parameter at coexistence remains constant independent 
of ionic strength, but in Figs. |^ and |l^, a weak de- 
pendance of the order parameter with ionic strength is 
seen. In Fig. |l^, the ionic strength dependance of the 
nematic order parameter at coexistence is plotted as de- 
duced from both x-ray diffraction and birefringence mea- 
surements. The change in ionic strength from 5 mM to 
110 mM corresponds to an L/D e s for the rods chang- 
ing from ss 40 to « 85. As the effective aspect ratio for 
our rods approaches the long rod limit, L/D c ff > 100, 
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FIG. 7: Concentration dependence of the nematic order pa- 
rameter, (a) is at 10 mM, and (b) is at 110 mM ionic strength 
and pH 8.2. Squares (■) are from the interparticle interfer- 
ence peak, and open circles (o) are results from the intraparti- 
cle peak. The solid lines shown are for a scaled particle theory 
for charged semi-flexible rods. 



the coexistence order parameter decreases, approaching 
the theoretically predicted value of S = 0.55, as calcu- 
lated by Chen for long semi-flexible rods with a length 
to persistence length ratio, L/p = 0.4 111). Even though 
the persistence length of fd virus is more than twice its 
contour length, and thus can be considered fairly rigid, 
all of our co-existing samples had a nematic order pa- 
rameter significantly lower than the Onsager prediction 
of S — 0.79 as measured by both diffraction and birefrin- 
gence. 

Stroobants et al., have shown that there is an ad- 
ditional electrostatic twisting factor which acts to mis- 
align adjacent particles and may effect the nematic or- 
der parameter at coexistence j45|. This effect scales as 
h = K~ 1 /D c g where k _1 is the Debye screening length. 
The effect of h on the coexistence concentrations of the 
system is predicted to be small jl5| , and the predicted de- 
crease in the nematic order parameter is only 2.4% when 



9 




10 15 20 25 

Concentration [mg/ml] 



FIG. 8: Concentration dependence of the order parameter 
of the nematic phase co-existing with the isotropic phase as 
determined from the intraparticle peak (open circles o) and 
the interparticle peak (solid squares ■). Increasing concen- 
tration of the coexistence concentration is due to increasing 
ionic strength ^ . The solid line is a linear fit to the combined 
sets of data and is presented as a guide to the eye. 
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FIG. 10: Ionic strength dependance of the order parameter 
of the nematic phase co-existing with the isotropic phase as 
calculated by x-ray diffraction measurements (solid squares 
■ ) and birefringence measurements (open triangles A). Solid 
line shows order parameter predicted by scaled particle theory 
for charged semi-flexible rods as in Fig. (?]. 
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FIG. 9: Comparison of measured birefringence An/c to the 
deduced x-ray order parameter S. Open shapes are from in- 
traparticle interference peak measurements. Closed shapes 
are from interparticle interference measurements. The equa- 
tion of the fitted line is An/c = (3.8 ± 0.3)S-(0.11 ± 0.19) 
where An/c is in units of 10 ml/mg. 



increasing ionic strength from 5mM to HOmM, whereas 
we measure an decrease in the order parameter of ~ 15%. 
Nevertheless, the observed trend of decreasing order pa- 
rameter with increasing ionic strength indicates that the 
electrostatic twisting effect might be influencing the sys- 
tem in addition to the changing effective aspect ratio of 
the charged virus. It is also important to note that be- 
low about 10 mM, the concentration of the virus may 
begin to have an effect on the ionic strength of the so- 
lution in that the concentration of the virus counterions 



may act to increase the ionic strength. Overall, we ex- 
pect and observe better agreement with the theoretical 
predictions for the nematic order parameter at high ionic 
strength. 



VI. CONCLUSION 

To summarize, we have observed, as predicted by Ham- 
ley, that the method of using x-ray diffraction to calcu- 
late the oricntational distribution function is insensitive 
to the details of the form of distribution function used. 
Nevertheless, we were able to rule out the function used 
by Oldenbourg et al. because we could qualitatively see 
that models created using this function did not fit the 
data at low concentration and at high angle \& from the 
equator on the detector. The Onsager and Gaussian trial 
angular distribution functions fit the angular distribu- 
tion of both the intraparticle and interparticle diffraction 
peaks equally well and returned similar values for the ne- 
matic order parameter. The concentration dependence 
of the nematic order parameter at high ionic strength, or 
large L/D c g, as determined from both the interparticle 
and intraparticle scatter agrees with that predicted by a 
scaled particle theory of charged semi-flexible rods. At 
low ionic strength, theoretical predictions qualitatively 
reproduce the concentration dependance of the order pa- 
rameter. Similar agreement of the concentration depen- 
dence of nematic ordering to Onsager's theory has been 
measured for the system of persistence lengthed DNA 
fragments using neutron scattering |l7], [24J]. This simi- 
larity demonstrates the universality of Onsager's theory 
and its applicability to semi- flexible systems. 
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The nematic order parameters derived from both in- 
terparticle and intraparticle scatter return similar results, 
implying that it is sufficient to use the easier, one dimen- 
sional analysis of the interparticle interference peak to 
calculate nematic order parameters as has been done for 
many years for thermotropic liquid crystals. It has also 
been shown that the relationship between the birefrin- 
gence and the nematic order parameter as calculated by 
x-ray diffraction is linear. From this relationship the sat- 
uration birefringence of fd was calculated. Subsequently, 
the order parameter can also be obtained simply by mea- 
suring the birefringence of a sample of nematic fd and 
rescaling it by the saturation birefringence. We note that 
the birefringence measurements were much less repeat- 
able than diffraction measurements, as can be observed 
by the large error bars throughout the entire range of 
data shown in Fig. ||. 

At high ionic strength, or large effective aspect ratio, 
we observed that the order parameter of the nematic 
phase coexisting with the isotropic phase was S » 0.6, 
close to the theoretically predicted value for semi-flexible 
rods and significantly lower than the theoretical value of 
S = 0.79 for rigid rods. With decreasing ionic strength 
however, a weak systematic increase in the nematic co- 
existence order parameter was found. This is consistent 
with both an decrease in the twist parameter and a de- 
viation of L/D e g from the long rod limit. In order to 
fully understand the interactions which are producing 
the nematic phase diagrams, particularly at lower ionic 
strength where L/ D c g is small new theories and simula- 
tions need to be developed which include a more complete 
picture of the complicated electrostatic interactions. 
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APPENDIX A: X-RAY DIFFRACTION 
ANGULAR ANALYSIS 

If we assume we are at sufficiently high scattering an- 
gle where intensity variations due to interparticle inter- 
actions are negligible then S(q r ,q z ) = 1, and we can 
measure the orientational distribution function from in- 
traparticle interference by comparing it to a simulated 
scatter created from the evaluation of Eq. ||. To evalu- 
ate Eq. H a three dimensional model for the single rod 
form factor I s (q r , q z ) was developed. 

A long rod Fourier transforms as a disk of thickness 
2n / L oriented perpendicular to the long axis of the rod. 



Because of the helical periodic structure along the long 
axis of fd, the Fourier transform of a single fd consists of a 
series of disks separated by a distance proportional to the 
reciprocal of the period jl6| . This is shown schematically 
in Fig. [ll]a. The radial intensity along these disks is a 
summation of Bessel functions whose exact form depends 
on the structure of the rod. When projected onto a screen 
these disks are visible as layer lines. The images shown 
in Figure show the zeroth and ± first layer lines. For 
our model, the radial intensities of the disks were approxi- 
mated by the scattered intensities along the middle of the 
zeroth and ± first layer lines, / m , of our most aligned ne- 
matic sample, times the radius q r at which that intensity 
is located and the width, a, of the Gaussian ODF, 

I s = I m \f2Tiaq r . (Al) 

\Z2naq r is the disorientation correction term. For a small 
amount of disorientation of rods, the radial intensity de- 
creases as l/q r . The effect of the disorientation is il- 
lustrated in Fig. [TT|b. This approximation method was 
developed by Holmes and Leigh, and is valid if the sam- 
ple from which the I m is taken was well aligned |37| . The 
nematic order parameter of our most aligned sample was 
S = 0.96 as measured from the interparticle interference 
peak. 

In order to model diffraction from a nematic phase of 
fd, this single particle scattered intensity is multiplied 
by a test ODF and integrated over all possible angles 
of orientation, as in Eq. |^. The intersection of the re- 
sulting three dimensional nematic form factor and the 
Ewald sphere is then "projected" onto a two dimensional 
"screen" and a final two dimensional image is created, 
as shown in Fig. |IT|a. The width of the Ewald sphere 
was chosen to equal the energy uncertainty of the ex- 
perimental beam AE/E = 5 ± 1 X 10~ 4 . The screen 
image is then convolved with the Gaussian point spread 
function of the experimental x-ray beam approximated 
as exp[— r 2 /2a 2 ] with a = 0.0063 A. A series of two di- 
mensional images were made for different orientational 
distribution functions with different amounts of disor- 
der, examples can be seen in Fig. |^b. All collected data 
within the range of q r = 0.19 — 0.33Ais fitted to these 
simulated images to find the ODF which minimized a 
computed chi-squared value shown in Eq. |j| This q- 
range encompassed the innermost peak on each of the 
three layer lines visible in the intraparticle interference 
scatter. 

To measure the orientational distribution function 
from the interparticle peak, the method of Oldenbourg 
et al. was used. Because we measure the angular spread 
of only one diffraction peak, Eq. ^ simplifies to a one 
dimensional integral at constant q r : 

/(*) = J <P {8) I s {oj) sin uduj (A2) 

where is the angular intensity distribution along an 
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FIG. 11: (a)Schematic of the Fourier space of a single rod tilted at a slight angle, (b) Schematic showing how the intensity 
along the center of the layer lines decreases as q r when there is a small amount of angular disorder. Three rod axes (vertical) 
are labeled 1,2,3 along with their corresponding contribution to layer lines 0,1,2 (horizontal) as shown, (c) Schematic showing 
effect of the thickness of the form factor disks on the scatter with changing w. The right hand image in (c) is an enlargement 
of the equatorial intersection of the Ewald sphere and I s . 



arc drawn at constant radius, "F is measured from the 
equator on the detector film, <E>(#) is the angular distri- 
bution function of the rods and uj is the angle between the 
rod and the incoming beam. 'F, 9 and ui are related by 
cos# = cos ^ sin uj. Even thought it was originally used 
for analyzing intraparticle scatter, this equation is iden- 
tical to that used for analyzing thermotropic interparti- 
cle scatter, except that Oldenbourg's method includes a 
term which accounts for the length of the rod by defining 
the single rod scattering as I s (lu) = 1/ sinw, for small 8, 
where lo is the angle between the rod axis and the x-ray 
beam as illustrated in Fig. [TTj g. This 1/sinw propor- 
tionality comes from the understanding that the Fourier 
transform of a rod of finite length is a ring with a fi- 
nite thickness, and as uj decreases 1/ sin a; increases and 
more of the disk intersects the Ewald sphere and is sub- 
sequently projected onto the detector screen. 



Analysis done on intcrparticlc interference from ther- 
motropic liquid crystals typically defines I s (u>) = 1 
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3q ] . It has been previously shown through calcu- 
lations that neglecting the angular width when calculat- 
ing the order parameter from interparticle interference 
scatter results in inaccurate values for the nematic or- 
der parameter for S > 0.8 But, in our analysis 
we observed that changing I s from 1 / sin u> to one in the 
interparticle interference scatter analysis did not have a 
significant effect on the calculated value of the nematic 
order parameter, nor did the % 2 values reveal any in- 
formation as to which I s better describes the data. We 
chose to include the effect of rod length in our interparti- 
cle scatter analysis to be consistent with our intraparticle 
scatter analysis, which requires a knowledge of the rod 
length. 



APPENDIX B: SCALED PARTICLE THEORY 

To compare the experimental results for the order pa- 
rameter to the theory, we use the scaled particle expres- 
sion for free energy of hard rods as was developed by 
Cotter and coworkers jl7|, The main advantage of 
the scaled particle theory is that it takes into account 
third and all higher virial coefficients in an approximate 
way and leads to very good agreement with simulation 
results for the I-N coexistence pS] ]. Therefore this the- 
ory should be more adequate at describing data at higher 
concentrations of rods. We also note that the expression 
for the free energy (Eq. Bl) reduces to Onsager's second 
virial approximation for very long rods (L/D — > oo). The 
free energy derived by Cotter is : 



F(6,</>,a) 
Nk b T 

+n 2 ((5,a) 



ln(0)+ln(l-0) + <T($(o!)) 

2 



1 



1 



n 3 («5,a) 



1 - 



where 4> is the volume fraction of rods 
v V V6 4 



(Bl) 



(B2) 



The coefficients Il 2 and n 3 are given by the following 
expressions: 



n 2 (M) = 3 + ^-^($(a)), 



(B3) 
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and parameter S is the overall length over diameter ratio 
of the spherocylinder given by 6 = (L + D)/D. The func- 
tions £(a) is the excluded volume interaction between two 
rods as derived by Onsager 



£(«) = 



212(a) 
sinh 2 (a) 



(B5) 



The expression that accounts for the rotational entropy 
of the rods and the entropy associated with the loss of 
configurations due to confinement of the bending modes 
of the semi-flexible rods in the nematic phase has been 
derived by extrapolating between the hard rod and the 
flexible chain limits {§ fo[ [Si). In this paper the expres- 
sion obtained by Hentschke is used for numerical calcu- 
lations 



a(a, —) = ln(a) - 1 + ne 01 H (a - 1) 

p 6p 
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use Onsager approximation for the orientational distri- 
bution function <E>(a) (Eq. lot) and minimize the scaled 



particle free energy in Eq. Bl with respect to the pa- 



After the expression for the free energy is obtained, we 



ramcter a to find the order parameter of the nematic 
phase at different rod concentrations. To find out the 
concentrations of rods in the coexisting isotropic and ne- 
matic phases we solve the conditions for the equality of 
the osmotic pressure and chemical potential. 

To take into account the fact the rods are charged, 
instead of using the hard core diameter D in our cal- 
culations we use an effective diameter D c g ||. Strictly 
speaking this re-scaling procedure by D e g is valid only 
for densities at which the system is described by the sec- 
ond virial approximation, therefore our theoretical pre- 
diction has an uncontrolled approximation. Despite this 
fact the agreement between the theory and the experi- 
ments is quite satisfactory. It is worth mentioning that 
there have been recent effort to extend the validity of the 
scaled particle theory to include repulsive interactions, 
however this theory was not included in our calculations 
|52l |53|. 
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